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Aanleiding  tot  het  onderzoek 

In  werkgebied  I  van  het  TNO-DO  onderzoekprogramma  worden  rekenmodellen  ontwikkeld 
voor  de  beschrijving  van  constructiegeluidoverdracht  van  machines  naar  de  scheepsconstructie 
via  verschillende  typ)en  fundaties.  Bij  een  studie  naar  de  zogenaamde  ’doos-achtige’  fundatie 
(TPD-HAG-RPT-950092)  Week  dat  onvoldoende  duidelijk  was  wat  de  invloed  is  van  de  vorm 
van  het  contactvlak  tussen  de  fundatie  en  de  scheepsconstructie  op  de  geluidoverdracht. 

Doel  van  het  onderzoek 

Het  afleiden  van  analytische  uitdrukkingen  die  beschrijven  wat  de  invloed  is  op  de 
geluidoverdracht  van  de  vorm  van  het  contactvlak  tussen  een  ’doos-achtige’  fundatie  en 
scheepsconstructie. 

Werkplan 

Het  uitvoeren  van  een  analytische  en  numerieke  modellering  van  het  reele  deel  van  de 
effectieve  admittantie  voor  verdeelde  aanstoting  van  een  plaatachtige  constructie  via  een 
rechthoekige  rand.  Onderzoeken  van  de  invloed  van  de  lengte-breedteverhouding  van  de 
rechthoek  op  de  admittantie  voor  een  uniforme  spanningsverdeUng  (’nulde  orde’)  en  een  lineair 
verlopende  spanningsverdeUng  (’eerste  orde’).  De  resultaten  zullen  worden  vergeleken  met 
eerder  gevonden  resultaten  voor  een  cirkelvormig  contactvlak  (TPD-HAG-RPT-950024). 

Conclusies 

Het  contactvlak  mag  als  een  punt  worden  beschreven  zolang  de  omtrek  van  het  contactvlak 
kleiner  is  dan  de  buiggolflengte  in  de  ontvangende  plaatconstructie.  Bij  grotere  contactvlakken 
treedt  uitmiddeling  op,  zodat  de  effectieve  geluidoverdracht  vermindert.  Deze  vermindering  kan 
worden  gekwantificeerd  met  behulp  van  de  in  dit  rapport  beschreven  schattingsprocedures. 
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Preface 

The  present  report  concerning  the  influence  of  the  interface  geometry  on  an  extended  version 
of  mobilities  constitutes  a  part  of  the  work  within  the  mid-term  defence  research 
(DO-MLTP).  It  is  the  third  part  of  a  study  of  the  mobility  concept  employed  in  cases  where 
the  excitation  is  distributed  over  some  closed  contour  interface. 

The  author  is  grateful  to  Dr.  C.  deJong  and  Mr.  F.  v.d.  Knaap,  TNO  Institute  of  Applied 
Physics,  for  pointing  out  the  discrepancy  between  theoretical  and  experimental  results  which 
initiated  the  work. 

The  research  has  been  sponsored  by  the  TNO  Defence  Research  Organization  under  the  mid¬ 
term  plan  (contract  A95/KM/133) 


Loughborough  University  of  Technology  in  April,  1996 

Bjom  A.T.  Petersson 
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1  Introduction 

In  some  recent  studies  [1],  [2],  the  dynamic  interaction  between  a  superstructure  and  the 
recipient  for  the  vibrational  energy  has  been  considered.  In  particular,  the  cases  with 
interfaces  which  are  large  in  comparison  with  the  governing  wave  length  of  the  recipient 
require  detailed  consideration  in  order  to  circumvent  the  often  rather  cumbersome  mixed 
boundary  value  problem  implied.  In  [3],  the  important  conclusion  was  drawn  that  for  a  large 
class  of  stmctural  configurations  pertinent  to  the  shipboard  structural  acoustics  transmission 
problem,  essentially  the  zero  and  first  order  stress  distributions  suffice  for  the  modelling.  It 
was  argued,  moreover,  that  for  not  too  extreme  geometrical  shapes  of  the  interface,  the 
annular  interface  constitutes  the  generic  case.  In  subsequent  experimental  work  [4]  it  has 
been  observed,  however,  that  deviations  from  calculated  results  were  obtained,  below  and  in 
the  transition  range  from  point-like  to  surface-like  excitation  for  a  rectangular  interface  and 
the  question  was  posed  whether  those  deviations  could  be  explained  by  a  shape  dependence. 
The  interface  problem  underlying  the  aforementioned  study,  of  course,  is  of  central 
importance  for  superstructures  such  as  machinery  seatings  which  are  installed  on  recipients 
for  the  vibratory  energy,  e.g.  the  bottom  structure  of  ships.  Interface  mobilities  of  plate-like 
structures  have  been  derived  [5]  for  the  case  of  annular  shapes  which  show  that  for  both  the 
zero  (transverse  force)  and  first  order  (lateral  force  or  moment)  excitation  distributions,  the 
real  parts  equal  those  found  for  ordinary  point  excitation  when  the  governing  wavelength  of 
the  recipient  is  larger  than  the  circumference  of  the  interface  and  decrease  markedly  at  high 
frequencies.  The  imaginary  part  of  the  zero  order  interface  mobility  is  small  and  negative  and 
can  be  termed  a  negative  spring  while  that  of  the  first  order  interface  mobility  is  stiffness 
controlled  in  the  region  of  point-like  behaviour.  Although  the  quantitative  details  may  vary 
with  the  shape  of  the  interface  the  overall  character  will  remain  the  same.  This  means  that 
with  a  mass  controlled  superstructure,  the  transmission  may  be  enhanced  due  to  matching  but 
no  system  resonance  will  occur  for  transverse  force  excitation.  If,  on  the  other  hand,  the 
dynamic  characteristics  of  the  superstructure  are  stiffness  controlled,  a  system  resonance  may 
occur  because  of  the  different  signs  of  the  imaginary  parts  involved.  Such  a  system 
resonance,  of  course,  can  be  very  pronounced  but  then,  affect  a  limited  frequency  region. 
For  a  lateral  force  or  moment  excitation,  a  mass-controlled  behaviour  of  the  superstructure 
will  unavoidable  lead  to  a  system  resonance  whereas  a  stiffness-controlled  behaviour,  at 
worst,  lead  to  a  broad  matching.  In  the  experimental  comparison  referred  to  above,  however, 
the  deviations  showed  a  reduced  power  transmission  which  can  be  explained  by  a  markedly 
altered  real  part  of  the  interface  mobility  in  the  rectangular  case. 

In  order  to  investigate  such  a  shape  dependence,  the  transmission  of  vibrational  power  at  a 
rectangular  interface  is  addressed.  The  theoretical  analysis  is  carried  out  for  both  a  primary 
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transverse  force  excitation  and  a  primary  lateral  force  or  moment  excitation.  As  in  the  case  of 
an  annular  interface,  a  plate-like  recipient  is  assumed. 
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2  Analysis 

2.1  Zero  order  excitation  distribution 

Consider  a  thin  homogeneous,  linearly  elastic  plate  subject  to  an  excitation  distributed  over  a 
strip  with  an  overall  rectangular  form  as  depicted  in  Figure  1.  Since  the  primary  objective  of 
the  present  study  is  to  clarify  the  influence  of  the  shape  of  the  interface  for  the  power 
transmission,  the  real  part  of  the  mobility  constitutes  the  appropriate  quantity.  Owing  to  the 
fact  that  the  dimensions  of  the  interface  can  be  larger  than  the  wavelength  of  the  free  flexural 
wave,  however,  the  mobility  concept  must  be  extended.  This  is  achieved  by  considering  the 
so-called  effective  interface  mobility,  defined  in  [3].  As  discussed  in  e.g.  [6],  a  proper 
analytical  model  in  the  shipboard  context,  is  realized  by  a  recipient  plate  of  infinite  extent. 


Figure  1 :  A  thin  plate  excited  over  a  rectangularly  shaped  strip. 

For  the  sake  of  physical  rigour,  the  strip  is  assumed  to  be  of  finite  width  and  the  stress  to  be 
uniformly  distributed  both  along  and  across  the  strip-like  interface,  see  Figure  2. 
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Figure  2:  Cross-section  of  strip  distributed  excitation. 

Therefore,  let  the  zero  order  stress  be  given  by, 

<T=^- 

AK{a  +  b)  (1) 

For  a  thin  infinite  plate,  the  velocity  due  to  a  distributed  excitation  can  be  written  [7] 
Y~\\a(x,.y,)E(k,R)dS 

s  .  (2) 

where  7°°  =  (oj SB'kg  is  the  ordinary  point  force  mobility  of  an  infinite  plate  and 

E{kR)  =  H^^\kR)-H^^\-ikR)  _  (3) 

Following  [3],  the  real  part  of  the  interface  mobility  can  be  established  by  considering  the 
far-field.  For  large  distances  between  the  excitation  and  an  arbitrary  response  position,  the 
propagation  function  above  simplifies  into 


EikR)  = 


2_-i{kR-K/4) 


7tkR 


(3a) 


Moreover,  under  such  conditions, 

R^  =  ~  2x^x^  +x^  +y^  -  2y^y^  +  y^ 

which  can  be  approximated  by 

R^  ~  /?o(l“2x^cos0//?o  “2y^sin0/i?o) 


(4) 
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R~  -;r^cos0-y^sin6 


This  means  that  since  essentially  the  phase  is  sensitive  to  the  small  geometrical  perturbations, 

.  f 

^  4A(a  +  Z7)'\jrol:i?o  [44 

a-A  b-A  1 


-  J  Je-'*(^o-^.cos0-y,sin0)^^^^^ 


-a+A-b+A 


The  integrations  are  readily  obtained  as 

00 I  ^  [r2sin(^fiicos0)1  r2sin(fci>sin0) 

Vf/vAjL/)  =  Y  -  I - >  -  - 

4A{a  +  b)  '1|;z^/?o[L  kcosO  J[  ^sin0 

r  2sin{/:(a- A)cos0)1  f  2sin(^(Z7- A)sin0) 


fccos0 


^sin0 


Since  sin(a  -  P)  =  sin  a  cos  /3  -  sin  /5  cos  a ,  Eq.  (6)  can  be  rewritten  as 


v(/^,0)  =  F“ 


Fgg  1  2  J  sin(^acos0)'|  rsin(fc^?sin6) 


A(a  +  Z?)  [L  kcosO  JL  fcsin0  J 

sin(^a  cos  0)  cos(^A  cos  0)  -  sin(M  cos  0)  cos{ka  cos  0) 
kcosO 

sin(^Z7sin0)cos(A:Asin  0)  -  sin(A:Asin  0)cos(tosin  0) 
fcsin0 


where,  for  small  A , 

00  2  r  sin{^acos0) 

v(Fo,0)  =  F  — -  -  czcos(^Z?sin0) — - - 

(a  +  b)  '\i7rfci?oL  kacosO 


+bcos{kacos9) 


sm{kbsm6) 
kb  sin  6 


Hence,  for  kA  less  than,  say,  0.1,  the  farfield  is  independent  of  the  width  of  the  strip 
interface. 
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A  power  balance  can  now  be  employed  to  determine  the  real  part  of  the  interface  mobility, 


W  =  Re[7o]  =  CBR^m" \\v{RQ,e)fde 


yielding 


^  ^  ff  m  ■  r^^sm{kacosd)  ,  ,sm(^^sin0) 

Re[}n]  = - T  \  a cos( kb sm 6) — - -  +  bcos(kacosd) -  dd 


2Ti{a  +  by 


kb  sin  G 


Let  p  =  4{a  +  b)  be  the  perimeter,  see  Figure  1,  then  the  alternative  form  of  the  interface 
mobility  is  obtained  as 


^  ;  p 


f  vSin  kp—cos6  .  ^  sin  kp—sinG 

Re[7o]  =  — 7”  f  —cos  kp—sinG  — — - - ^  +— cos  kp—cosG  - -r^ - 

^  0  ^  ^  ^  ^  kp—cosG  P  \  P  J  kp—sinG 

P  P 


To  get  an  analytical  overview,  the  integrands  can  be  expanded  in  terms  of  sin  G  and  cos  G 
which  for  ka  «  1  and  kb  «  1  leads  to, 

■^oo  2/r 

ReTO-.  ,  ^,,2  i {(“+>’? 

l7Z{a  +  b) 


7  2  b^  2  1 

-2{a  +  b)k^  — (3^  +  ai)cos^0  +  — (3a  +  Z7)sin^0  >dG 
6  6 


and  as  a  limit  one  finds, 

"*2  ”  2 
ReKJ^y"  l-L(a  +  i,)2  =y-  l-L[4[a  +  b)f 
6  96 


This  is  to  be  compared  with  that  of  the  annular  case  [5]  i.e.. 


Re[7ol=r-(y„()cro)f  l-^(2»of  =r-  l-^{2»b) 

OTT  oU 


TNO-report 


TPD-HAG-RPT-960094 


April  1996 


page  1 1 


and  it  is  seen  that  the  decrease  of  the  interface  mobility  compared  with  the  ordinary  point 
force  mobility  is  slightly  larger  in  the  annular  case  than  in  the  case  of  a  rectangular  strip  for 
the  same  perimeter  as  schematically  illustrated  in  Figure  3. 


Helmholtz  number,  kp 


Figure  3:  Schematic  comparison  of  asymptotic  behaviour  of  annular  ( - )  and 

rectangular  ( - )  interface  mobility  for  the  same  perimeter  p. 


Rewriting  the  mobility  as  a  function  of  Helmholtz  number,  ka,  means  that  the  asymptotic 
behaviour  for  small  Helmholtz  numbers  is  captured  by 


Ro[Yq]^Y- 


1- 


(11) 


which  can  be  used  in  a  comparison  of  the  influence  of  different  shapes  of  the  rectangular 
strip. 

In  Figure  4  is  displayed  the  interface  mobility  versus  Helmholtz  number  for  various  aspect 
ratios  b/a.  It  is  clearly  seen  that  the  onset  of  the  reduction  is  controlled  by  the  aspect  ratio 
and  that  a  decade  increase  in  the  aspect  ratio  means  a  shift  of  a  decade  downwards  in 
Helmholtz  number  of  the  onset. 
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Helmholtz  number,  ka 

Figure  4:  Interface  mobility  versus  Helmholtz  number  based  on  side  length  for  some  aspects 

ratios:  bla  —  \  (  ),  bla-l  ( . ),  bla  =  1  ( - ),  bla  =  5  ( - )  and 

b!  a  =  \Q  Slope  of  characteristic  beam  mobility  ( - ) 

In  order  to  assess  the  applicability  of  the  approximation  in  Eq.  (1 1),  the  results  are  compared 
with  those  obtained  numerically  using  Eq.  (8).  In  Figure  5  is  shown  the  transition  region  and 


Helmholtz  number,  ka 


Figure  5:  Asymptotic  behaviour  of  the  interface  mobility  for  some  aspects  ratios:  b/ a  =  l 
( - ),  b/a  =  2  ( . ),  bla  =  3{ - ),  b/a  =  5  ( - )  and  b/a  =  lO 
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it  emerges  that  the  approximation  for  small  Helmholtz  numbers  is  applicable  in  the  range 
below  ka  <  1/(1  +  bl  a)  but  deviates  markedly  from  the  integrated  results  above  this  limit. 

A  perhaps  more  adapted  way  of  displaying  the  variation  in  interface  mobility  with  aspect  ratio 
is  as  function  of  Helmholtz  number  based  on  the  perimeter,  used  in  Figure  6. 


Figure  6:  Interface  mobility  versus  Helmholtz  number  based  on  perimeter  for  some  aspects 

ratios:  bla  =  \  ( - ),  b/a  =  2  ( . ),  bla  =  3  (-  -  -  -  -),  bla  =  5  ( - )  and 

b/a  =  10  ( - ). 

Here,  as  indicated  in  Figure  3,  it  is  evident  that,  when  the  perimeter  equals  a  wavelength  of 
the  governing  wave,  the  spatial  matching  of  the  free  wave  and  the  excitation  distribution 
starts  to  deteriorate  and  the  real  part  of  the  interface  mobility  diminishes  irrespective  of  the 
aspect  ratio.  Hence,  the  tentative  conclusion  drawn  in  [5]  is  strongly  corroborated  for  this 
excitation  distribution  and  accordingly  the  perimeter-wavelength  criterion  appears  applicable 
for  the  onset  of  the  influence  of  the  distributed  excitation. 

The  more  or  less  deep  troughs,  of  course,  stem  from  the  fact  that  power  circulation  occurs 
between  the  sides  of  the  rectangular  interface.  Only  in  the  case  of  a  square-shaped  strip 
interface,  the  pattern  is  regular. 

Since,  as  illustrated  in  Figure  4,  the  envelope  for  the  upper  range  of  Helmholtz  numbers 
follows  a  square  root  law  it  is  possible  to  establish  a  simple  estimation  procedure  as  depicted 
in  Figure  7.  For  the  upper  range,  the  estimates  will  be  conservative.  By  shifting  the  onset  of 
the  inverse  square  root  dependence  to  kpjlK  =  1/2 ,  the  estimation  in  the  upper  range  is 
somewhat  improved. 
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Figure  7:  Estimation  procedure  for  the  zero  order  real  part  of  interface  mobility  of  closed 
contour  strip  interfaces. 
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2.2  First  order  excitation  distribution 

For  this  case  as  well,  the  strip  is  assumed  to  be  of  finite  width  but  the  stress  to  be  linearly 
distributed  in  the  jc-  or  y-directions  and  to  have  a  uniform  distribution  along  the  remaining 
strips,  see  Figure  8. 


(7 


Figure  8;  Cross-section  of  strip  distributed  moment  excitation. 


Therefore,  let  the  stress  take  the  form 

fsign(a);  -b<y<b,  \x\-a 


a{x,y)  = 


Mn 


4a(a/3rt-Z?)A 


;  -  a<  x<  a  ,  |y|  =  ^ 


(12) 


As  in  the  previous  sub-section,  a  power  balance  is  employed  and  leads  to, 

2n 

,MJ^  RelFI  =  CuRr^m 
2' 


-  Ik 

W  =  -\m/  Re[7]  =  CsRom"  l\v(Ro,Bfde 


0  ,  (13) 

from  which  the  real  part  of  the  first  order  interface  mobility  is  formally  found  to  obey, 

,2;r|  a  b 


Re[yi]  =  1 1  j  1  cos0-y, 


0  \-a—b 


de 


(14) 


The  inner  integrations  are  readily  obtained  as. 
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7=4/1 


bsm[kacos0  ) 


sm[kbsm6  ) 
kbsinO 


+  aco^kbsmB  ) 


sin[kacos0  )  cos{kacos0  ) 


(^acos0  kacosO 


.(15) 


The  real  part  of  the  first  order  interface  mobility  thus  becomes, 

ttOO  IKT"  1 


Rom 


27i{a/3  +  bY 


0  o 


— sin  (to  cos  0) 
a 


sin(tosin0) 
to  sin  6 


+cos(tosin0) 


sin(tocos0)  cos(tocos0) 
(to  cos  to  cos  0 


de 


,(16) 


which  also  is  independent  of  the  width  of  the  strip  interface  for  small  enough  kA  ■ 

As  for  the  uniform  stress,  the  integrand  can  be  expanded  in  terms  of  sin  Q  and  cos  6  to  get  an 
analytical  overview  which  for  to  «  1  and  kb  «  1  leads  to. 


Re[Fi] 


Y°°k^ 


2K(a/3  +  bY 


2k 

cos"  0 


——cos 0(0  +  al3)  b cos 0^0^  sin^  0  +  a~  cos^  0 j 


3 

-a  cos  01 


— cos^0  +  0^  sin^0 

v5 


>de 


and  after  integration,  the  asymptotic  behaviour  is  found  to  be 

2  r  2 

RerF ]  -4  F~  —  i  1 - — — (8  +  a / 0)  +  50^ (4  +  30 / a)  +  500a 

*  2[  18O(a/3  +  0)2‘-  ■ 

As  expected,  the  real  part  of  the  first  order  interface  mobility  converges  towards  that  of  the 
ordinary  point  moment  mobility  as  the  wavenumber  and  the  dimensions  decrease.  It  is  clear 
from  Eq.(17)  however,  that  the  dependence  on  the  aspect  ratio  is  more  involved  than  in  the 
case  of  a  uniform  stress.  Indeed,  the  asymptote  demonstrates  an  area  dependence.  The  real 
part  of  the  first-order  interface  mobility  for  the  rectangular  case  can  be  compared  with  that  of 
the  annular  case,  i.e.. 


TNO-report 


TPD-HAG-RPT-960094 


April  1996 


page  17 


Re[Fi]=  27° 

and  so  by  rewriting  Eq.  (17)  as, 

f,  (kafb/a 


riiCA'b)) 

1^  J 

k^ 

1  'O  J 

'  2 

4 

2 

158 

Re[7i]->7 


1- 


180(l/3  +  ^/fl)- 


-  [3(8  +  alb)  +  5ib/  af  {4  +  3b/a)  +  50b  /  a] 


(18) 


,(17a) 


it  is  seen  that  for  a=b  and  letting  a  correspond  to  the  radius  /q  ,  the  rectangular  case  implies 
a  smaller  mobility; 


Re[7i]->^7“'|1 


In  contrast  to  the  case  with  a  uniform  stress  excitation,  it  is  seen  that  the  decrease  of  the 
present  interface  mobility,  normalised  with  respect  to  the  ordinary  point  moment  mobility,  is 
slightly  smaller  in  the  annular  case  than  in  the  case  of  a  rectangular  strip  for  the  same 
perimeter.  This  schematically  illustrated  in  Figure  9. 


Helmholtz  number,  kp 


Figure  9:  Schematic  comparison  of  asymptotic  behaviour  of  annular  (■ 
shaped  ( - )  first  order  interface  mobility  for  the  same  perimeter  p. 


■)  and  square- 
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In  Figure  10  is  displayed  the  interface  mobility  versus  Helmholtz  number,  based  on  the  short 
side  length,  for  various  aspect  ratios  b/a.  It  is  seen  that  the  onset  of  the  reduction  is 
controlled  by  the  aspect  ratio  and  that  a  decade  increase  in  the  aspect  ratio  means  a  shift  of  a 
decade  downwards  in  Helmholtz  number  of  the  onset.  One  can  note,  however,  that  the 
signature  in  the  upper  range  is  more  regular  than  that  of  the  uniform  stress. 


Helmholtz  number,  ka 

Figure  10:  Interface  mobility  versus  Helmholtz  number  based  on  side  length  for  some 

aspects  ratios:  bla  =  l  ( - ),  b/a  =  2  ( . ),  b/a  =  3  { - ),  b! a  - 5  ( - ) 

and  ij/a  =  10  ( - ). 

In  Figure  1 1  is  visualised  the  range  of  validity  of  the  asymptotes  developed  for  small 
Helmholtz  numbers.  Again,  the  approximation  given  by  Eq.  (17a)  is  applicable  in  the  range 
below  ka  <\l{\  +  bl  a)  and  deviates  markedly  from  the  integrated  results  in  the  range  above. 
The  adapted  presentation  of  the  first  order  interface  mobility  as  function  of  Helmholtz 
number  based  on  the  perimeter,  is  given  in  Figure  12.  It  is  again  observed  that  the  onset  of 
the  reduction  in  real  part  occurs  when  the  perimeter  equals  the  wavelength  of  the  governing 
wave. 

Owing  to  the  fact  that  the  first  order  excitation  distribution  is  only  single  symmetric,  it  is 
necessary  to  numerically  analyse  also  the  influence  of  the  shape  of  the  interface  for  aspect 
ratios  smaller  than  unity.  In  Figure  13  are  shown  the  real  parts  of  the  moment  interface 
mobilities  for  some  aspect  ratios  smaller  than  unity. 
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Helmholtz  number,  ka 

Figure  11:  Asymptotic  behaviour  of  the  interface  mobility  for  some  aspects  ratios:  b! a  =  \ 
( - ),bla  =  2  ( . ),  b/a  =  3( - ),bla-5{ - )andb/a  =  lO  ( - ). 


1 


Figure  12:  Interface  mobility  versus  Helmholtz  number  based  on  perimeter  for  some  aspects 

ratios:  bla  =  \  ( - ),  bfa  =  2  ( . ),  bla-3  (-  - - ),  b/a  =  5  ( - )  and 

b/a  =  lO 
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Helmholtz  number,  kb 

Figure  13;  Interface  mobility  versus  Helmholtz  number  based  on  perimeter  for  some  aspects 

ratios:  b/a  =  l  ( - ),  bla=l/2  ( . ),  b/a  =  1/3  ( . ),  b/a  =  1/5  ( - )  and 

b/a=  1/10  ( - ). 

It  is  seen  that  when  the  aspect  ratio  becomes  smaller  than  unity,  i.e.  the  interface  becomes 
long  in  the  direction  perpendicular  to  the  moment,  the  pattern  of  constructive  and  destructive 
is  irregular,  cf.  Figure  10.  Although  such  a  trend  is  indicated  also  for  the  zero  order 
distribution  it  is  comparatively  pronounced  for  the  first  order  distribution  owing  to  the 
basically  elliptical  field  in  the  plate.  The  irregularity  thus  stems  from  the  fact  that  here, 
sections  of  the  interface  stress  can  interfere  when  the  long  direction  of  the  interface  is 
perpendicular  to  the  major  axis  of  the  elliptical  field. 

The  asymptotic  behaviour  for  small  Helmholtz  numbers,  however,  is  still  valid  as  seen  in 
Figure  14.  Of  course  it  must  be  noted  that  when  Helmholtz  number  is  based  on  half  the  side 
length,  b,  the  region  of  validity  is  given  by  kb  <\l{\  +  a/ b)- 
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Helmholtz  number,  kb 


Figure  14:  Asymptotic  behaviour  of  the  interface  mobility  for  some  aspects  ratios:  b/  a  =  l 
( - ),  bla  =  1/2  ( . ),  b/a=ll3  { . ),  b/a  =  1/5  ( - )  and  b/a  =  1/10 

For  the  sake  of  completeness,  the  interface  mobility  is  displayed  also  as  function  of 
Helmholtz  number  based  on  the  perimeter,  see  Figure  15. 


Figure  15:  Interface  mobility  versus  Helmholtz  number  based  on  perimeter  for  some  aspects 

ratios:  b/a  =  1  ( - ),  b/a^  1/2  ( . ),  b/a=  1/3  ( - ),  b/a=  1/5  ( - )  and 

b/a=l/10 


TNO-report 


TPD-HAG-RPT-960094 


April  1996 


page  22 


Accordingly,  when  the  aspect  ratio  falls  below  unity,  the  normalised  interface  mobilities 
versus  Helmholtz  number  based  on  the  perimeter  cluster  and  the  aspect  ratio  i.e.  the  shape 
becomes  more  or  less  irrelevant  provided  the  interface  is  long  in  the  direction  perpendicular 
to  the  moment  vector. 

As  for  the  zero  order  interface  mobility,  one  may  rather  easily  identify  the  overall  trends  of 
the  first  order  interface  mobility.  If  one  lets  the  length  of  the  rectangular  interface  become 
very  large  in  the  direction  of  the  moment,  then  the  problem  approach  a  two-dimensional 
counterpart  i.e.,  a  line  moment  on  a  plate.  In  such  a  case  the  slope  is  proportional  to  the 
square  root  of  frequency.  As  the  wavelength  approach  the  length  of  the  short  side  of  the 
rectangular  interface,  the  spatial  mismatch  between  the  excitation  and  the  wave  field  increases 
markedly  and  a  reduction  in  the  power  transmission  is  obtained.  In  this  region  therefore  the 
slope  becomes  inversely  proportional  to  the  square  root  of  frequency,  as  demonstrated  in  [5]. 
The  added  complication  with  an  estimation  procedure  for  the  first  order  interface  mobility 
hence,  is  the  transition  point  between  the  two  upper  regions.  As  discussed  in  conjunction 
with  Eq.  (17),  a  direct  relationship  does  not  exist  between  the  mobility  and  the  aspect  ratio 
nor  the  perimeter  but  a  combined  dependence  on  aspect  ratio  and  the  area  enclosed  is  found. 
By  focusing  on  the  range  of  the  aspect  ratio  bla<5,  however,  a  reasonably  straight¬ 
forward  procedure  can  be  established  from  the  critical  point  where  the  length  of  the  short  side 
of  the  rectangular  interface  equals  a  wavelength.  This  means  that  the  estimation  procedure 
can  be  summarised  as, 

kpf27:<l\  b/a<5 

2p 

Y°°-^kp  ;  1  <  kpl2K  <Kbl a\  1  <  Z?/a  <  5 
P 

3  1 

;  7cbla<kpl2K\  b/a<5 

P  (19) 

In  Figure  16,  the  results  using  the  estimation  procedure  are  compared  with  the  calculated 
results  for  the  limiting  case  b/ a  =  5.  It  is  seen  that  the  overall  signature  of  the  first  order 
interface  mobility  is  captured  by  the  simple  procedure  although  minor  under-  and 
overestimations  are  obtained  in  limited  regions. 
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Helmholtz  number,  kp/2Tt 


Figure  16:  Estimation  procedure  for  the  real  part  of  first  order  interface  mobility  of  closed 

contour  strip  interfaces,  displayed  for  b/a  =  5.  ( - )  Re[Yj(6/a=5)],  ( . )  generating 

curves. 

This  means  that  when  the  aspect  ratio  is  in  the  range  of  1  <  <  5,  the  first  order  interface 

mobility  has  three  distinct  regions,  a  first  in  which  the  ordinary  point  moment  mobility  gives 
an  adequate  description,  an  intermediate  in  which  the  principal  behaviour  is  that  of  a  moment 
excited  beam  and  finally,  an  upper  region  in  which  the  interface  mobility  is  markedly  reduced 
compared  with  plate  moment  mobility.  As  the  aspect  ratio  falls  below  unity,  the  intermediate 
region  vanishes  and  the  first,  constant  region  joins  the  third,  steep  region.  For  aspect  ratios 
larger  than  five  the  uppermost,  steep  region  can  be  omitted  if  the  range  of  interest  is  below 


the  upper  transition  point,  approximately  given  by  kp  ~  2-^n^b  /  a  ■ 
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3  Discussion  and  concluding  remarks 

The  influence  of  interface  geometry  on  the  transmission  of  vibrational  power  from  a 
superstructure  to  a  recipient  has  been  studied.  The  limited  effect,  previously  heuristically 
inferred,  is  theoretically  confirmed.  It  can  thus  be  concluded  that  as  long  as  the  perimeter  is 
less  than  the  wavelength  of  the  governing  wave  in  the  recipient,  the  ordinary  point  quantities 
are  valid.  For  larger  Helmholtz  numbers  where  the  perimeter  is  larger  than  the  wavelength, 
the  real  part  of  the  interface  mobility  exhibits  more  or  less  pronounced,  rounded  maxima  and 
sharp  minima  set  by  the  constructive  and  destructive  interferences  at  the  enclosed  area  of  the 
plate.  The  envelopes  to  the  zero  and  first  order  interface  mobilities,  in  the  case  of  a  hollow 
rectangular  interface,  are  those  obtained  for  an  annular  one  [5]. 

In  this  study,  focusing  on  the  vibration  transmission,  the  imaginary  part  of  the  interface 
mobility  is  not  treated  explicitly.  This  is  so,  foremost,  because  at  large  Helmholtz  numbers, 
this  part  is  sensitive  to  the  details  of  the  interface  geometry  as  well  as  the  stress  distribution 
cf.  [8]  and  thus  any  result  is  valid  merely  for  a  specific  configuration.  Qualitatively, 
however,  the  behaviour  for  large  Helmholtz  numbers  are  those  of  the  annular  case  when  the 
aspect  ratio  of  the  rectangular  interface  is  not  extreme,  typically  less  than  3,  and  tends  to  that 
of  the  strip  mobility  [9]  for  large  aspect  ratios.  This  means  that  in  a  broad-band  sense,  the 
transmission  is  controlled  by  the  real  part  of  the  interface  mobility.  Therefore,  it  is  not  likely 
that  the  overall  deviations  from  calculated  results  noted  in  the  experimental  work  on  a  large 
freely  suspended  plate  are  due  to  an  interface  shape  dependence  [4].  The  experimental  results 
reported  therein  for  a  rectangularly  shaped  box  with  an  aspect  ratio  of  1.5,  however,  showed 
no  pronounced  minima,  as  is  predicted  for  an  annular  or  square-shaped  interface.  As 
displayed  in  Figure  6,  this  feature  is  reduced  for  aspect  ratios  larger  than  unity  which 
demonstrates  that  such  destructive  interference  effects  are  indeed  dependent  upon  the  shape 
of  the  interface. 

The  zero  order  interface  mobility  is,  qualitatively  as  well  as  quantitatively,  practically 
independent  of  shape  provided  the  perimeter  is  retained  and  the  aspect  ratio  is  close  to  unity. 
As  the  aspect  ratio  grows  large,  the  interface  mobility  approaches  that  of  a  point  excited 
beam,  again  largely  irrespective  of  shape.  Accordingly,  these  two  trends  constitute  a 
conservative  basis  for  an  estimation  procedure. 

In  contrast,  the  first  order  mobility  for  a  rectangular  interface  must  be  treated  differently  for 
aspect  ratios  larger  and  smaller  than  unity  respectively.  In  the  former  case,  the  interface 
mobility,  at  large  Helmholtz  numbers,  tends  to  the  corresponding  beam  mobility  (moment 
excitation)  whereas  in  the  latter,  the  'moment'  interface  mobility  decreases  as  the  force 
interface  mobility,  cf.  [5].  This  means  that  two  different  bases  are  required  for  an  estimation 
procedure  with  respect  to  the  first  order  interface  mobility. 
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For  a  first  order  excitation  it  is  not  necessarily  beneficial  to  design  the  interface  long  in  the 
direction  of  the  moment  vector  whereas  it  is  always  advantageous  to  increase  the  length  in  the 
direction  across  the  moment  vector. 

For  both  the  zero  and  first  order  distributions,  there  are  more  or  less  pronounced  minima.  In 
principle,  such  minima  can  be  used  for  detuning  the  transmission  for  pronounced  tonal 
components  in  the  source  descriptor.  Although  only  the  case  with  an  aspect  ratio  of  unity 
results  in  a  regular  pattern  for  all  the  minima,  the  pattern  is  well  organised  also  for  large 
aspects  ratios  already  at  the  third  minima.  The  exception  is  the  first  order  interface  mobility 
where  a  long  interface  across  the  moment  vector  reduces  the  significance  of  the  local  minima 
and  establishes  an  irregular  signature. 
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Appendix  I: 

Symbols  and  notation 

B’ 

flexural  stiffness  per  unit  length 

E 

propagation  function 

F 

force 

H 

Hankel  function 

I 

integral 

J 

Bessel  function 

M 

moment 

R 

radius 

S 

area 

W 

power 

Y 

mobility 

a,  b 

half  side  length  of  interface 

c 

phase  speed 

i 

imaginary  unit 

k 

wavenumber 

m” 

mass  per  unit  area 

P 

perimeter 

V 

translational  velocity 

w 

rotational  velocity 

X.  y,  z 

Cartesian  co-ordinates 

A 

width  of  strip  interface 

a,l3 

auxihary  variables 

d 

polar  co-ordinate 

CO 

angular  frequency 

G 

normal  stress 

indices; 

B 

flexural 

ann 

annular 

rect 

rectangular 

r 

receiver 

s 

source 
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radius  from  interface  centre, 
net  external,  zero  order 
first  order 
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Geluidarm  Construeren,  ir.  J.C.  Vellekoop*. 

10.  TNO-TPD-TU  Delft,  Hoofdafdeling  Geluid,  Afdeling  Geluidarm 

Construeren,  dr.ir.  C.A.F.  de  Jong*. 

11.  TNO-TPD-TU  Delft,  Hoofdafdeling  Geluid,  Afdeling  Geluidarm 

Construeren,  dr.  B.A.T  Petersson*. 

12.  TNO-TPD-TU  Delft,  Hoofdafdeling  Geluid,  Afdeling  Geluidarm 

Construeren,  ing.  F.G.P.  van  der  Knaap*. 

13.  TNO-TPD-TU  Delft,  Hoofdafdeling  Geluid,  Afdeling  Geluidarm 

Construeren,  prof.dr.ir.  J.W.  Verheij*. 

14.  TNO-DO*. 

15.  TNO-TPD-TU  Delft,  Hoofdafdeling  Geluid,  Archief. 

16.  (t/m  18.)  Bibliotheek  KMA*** 

Indien  binnen  de  krijgsmacht  extra  exemplaren  van  dit  rapport  worden  gewenst  door  personen 

of  instanties  die  niet  op  de  verzendlijst  voorkomen,  dan  dienen  deze  aangevraagd  te  worden  bij 

het  betreffende  Hoofd  Wetenschappelijk  Onderzoek  of  indien  het  een  K-opdracht  betreft,  bij  de 

Directeur  Wetenschappelijk  Onderzoek  en  Ontwikkeling. 

*  Alijd  een  volledig  rapport. 

**  Indien  opdrachtgever  een  volledig  rapport,  anders  een  beperkt  rapport. 

***  3  volledige  rapporten  indien  het  rapport  ongerubriceerd  of  ongemerkt  is,  anders  1 

volledig  rapport. 


